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Abstract 

Anisotropic Alfvenic fluctuations with fc||/fc_L <^ 1 remain at frequencies much smaller than 
the ion cyclotron frequency in the presence of a strong background magnetic field. Based on 
the simplest truncation of the electromagnetic gyrofluid equations in a homogeneous plasma, a 
model for the energy cascade produced by Alfvenic turbulence is constructed, which smoothly 
connect the large magnetohydrodynamics (MHD) scales and the small "kinetic" scales. Scaling 
relations are obtained for the electromagnetic fluctuations, as a function of k± and /cy. Moreover, a 
particular attention is paid to the spectral structure of the parallel electric field which is produced 
by Alfvenic turbulence. The reason is the potential implication of this parallel electric field in 
turbulent acceleration and transport of particles. For electromagnetic turbulence, this issue was 
raised some time ago in [A. Hasegawa, K. Mima, J. Geophys. Res. 83 1117 (1978)]. 
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I. INTRODUCTION 



At large wave-numbers in the direction perpendicular to the ambient magnetic field, 
Alfven waves produce a significant compression of the plasma which results in the creation 
of a parallel electric field via the thermo-electric effect. The parallel electric field associated 
with a spectrum of kinetic Alfven waves leads to wave-particle interactions that can cause 
turbulent transport phenomena. This was pointed out already some time ago in Reference 
[1], in the context of the magnetospheric plasma. Starting from the drift- kinetics for the 
electrons, the authors of Ref.[l] derived a quasi-linear diffusion equation in both momentum 
and coordinate space, a so-called double-diffusion, showing the pivotal role played by the 
parallel electric field spectrum. In fact, their result is valid for any kind of low-frequency 
electromagnetic turbulence, but they used the properties of the kinetic Alfven wave to give 
order of magnitude estimates of various transport coefficients (density, momentum) and the 
electrons heating rate. The solar wind is a natural laboratory plasma where the spectral 
properties of electromagnetic turbulence can be probed with good accuracy down to scales 
below the ion Larmor radius. Turbulent Alfvenic fiuctuations observed in the solar wind 
have a solar origin. Hence, beside their intrinsic interests, solar wind measurements can also 
be used as a diagnostic tool to infer the properties of the turbulence closer to the sun, in 
the corona, where the plasma conditions are however different. For low- frequency Alfvenic 
fluctuations, the spectral energy density of the parallel electric fleld remains a small fraction 
of the electro-magnetic energy density. Hence, it is hardly a measurable quantity. However, 
since the parallel electric field is the actual force which mediates the wave-particle interaction, 
it is interesting to investigate its spectral structure in this range of spatio-temporal scales. 
This task is motivated by its potential relevance for turbulent acceleration and transport and 
is facilitated by more recent advances in the understanding of Alfvenic turbulence, mainly 
from theory and observations of the solar wind. 

The starting point is a minimal fiuid model of the non-linear dynamics of low-frequency 
kinetic Alfven waves, also capable of capturing the effect associated with the thermal Lar- 
mor radius of the ions, pi — vrji^ci with — ^jTi/rrii and uici — eBo/miC. This model 
can be obtained as a truncation of the electromagnetic gyrofiuid equations in a constant 
background magnetic field, by assuming a constant temperature for both the ions and the 
electrons, a constant background ion density and discarding all but the lowest two parallel 
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moments of the electron kinetics (see [2-5] and references therein). These gyrofluid equations 
are moments of the gyrokinetic equations. Gyrokinetics owe its name from being based on 
an averaging of the kinetic and Maxwell equations over the gyromotion of the particles and, 
hence, it applies in the limit of frequencies small compared to the ion cyclotron frequency, 
u <C cUci- Extensive reviews on gyrokinetics can be found in Refs.[6-8]. Gyrokinetics and 
their fluid moments are valid in the limit of small Larmor radius, pi <^ L, L being the 
macroscopic length scale of the plasma, and small anisotropic fluctuations, in the sense that 
k\\/k± ~ 6B/Bn -C 1, -Bo being the magnitude of the background magnetic field. As for 
reduced-MHD |9l4l l| . these ordering imply pressure balance in the direction perpendicular to 
the background magnetic field, such that the fast mode is ordered out. Moreover, because 
the gyroaverage procedure eliminates the cyclotron resonance, the only type of wave-particle 
interaction that remains possible is through the Landau resonance between the particles and 
the parallel electric force or the magnetic mirror one. On other central assumption of gyroki- 
netics, and hence of its fluid counterpart, is the small deviation of the distribution functions 
of the particles from a background distribution which is here a Maxwellian distribution. 

The electromagnetic gyrofluid equations are presented first in Section II. Then they are 
used, in Section III, to construct a model of the energy cascade of Alfvenic turbulence from 
the large MHD scales down to the small kinetic scales. The spectral structure of the parallel 
electric field produced by Alfvenic turbulence is studied in Section IV. Finally, possible 
implications of the existence of this parallel electric field are discussed in Section V. 

II. THE ELECTROMAGNETIC GYROFLUID MODEL 

The gyrofluid model considered in this work involves the electron density rif., the magnetic 
flux function = —Az, where A is the vector potential, and the electrostatic potential </>, z 
being the coordinate along the background field Boj2 5|. Adopting an MHD normalization 
(see the Appendix), these equations are 

atne+[0,ne]-V||J = O, (1) 

+ V||(pX-0) = 0. (2) 

The Poisson bracket is defined as [f,g] = z.V/ x S/g, where z is the unit vector along Bq, 
and V||/ = dzf + [ip, f] for any fields / and g. The quantity J = V^ip is the parallel current 



density and Ps — Csjojf^L is the normalized ion sound Larmor radius with Cg — ^jTjrni. 
This system is closed by the gyrokinetic Poisson equation, 

(ro - 1) , ... 

He = 2 ^' (3) 

Pi 

where Fq is an integral operator which describes the average of the electrostatic potential 
over a ring of Larmor radius pj. In Fourier space, Fo(6) simply becomes 

Fo(6) = e-%(6), (4) 

where h — pfkj^ and Iq is the modified Bessel function of the first kind. The above equations 
conserve the total energy given by 

E^l d^iiV^^f + plnl + 0(1 - Fo)0/p^]. (5) 

The model is valid for v\^lv\ <C 1, and moreover, since the effects associated with the 
electron inertia are not taken into account, it means that the conditions v\Jv\ ^ 1 and 
kx_de <^ 1 must also be satisfied, va being the Alfven speed and de the electron skin depth. In 
other words, the present model only describes the "kinetic" regime of the dispersive Alfven 
wave, not its "inertial" regime. 

The following comments are due. When a Fade approximant of the operator Fo(6) is 
used, which gives : 

Toib) - 1 ^ (6) 
the Foisson equation is converted into: 

(1 - p^VDrie = Vic/). (7) 

The model (l)-(3) with the Fade approximation (7), also including electron inertial effects, 
was studied in the context of collisionless magnetic reconnection in Ref.[5]. A well known 
effect of the electron inertia is the breaking of the frozen-in condition for the magnetic field. 

When k±pi <^ 1, relation (7) expresses the fact that the density is also equal to the 
plasma vorticity Ue — V±4>, in which case Eqs.(l)-(2) become equivalent to a low-/3 Hall- 
MHD model (see Appendix and e.g. Ref.[12]). Reduced Hall-MHD is an extension of the 
standard reduced MHD with the Hall and the electron pressure effect accounted for in the 
Ohms's law. The parallel component of the latter, i.e. Eq.(2), can also be written as 

^11 = -p'Vyne, (8) 



which shows the existence of a parallel electric field that arises from electron pressure gradient 
along the field lines. Electrons are here assumed to be isothermal. On the contrary, = 
is the standard reduced MHD Ohm's law. Equation ([8j) is a relation between the parallel 
electric field and the compressive density fiuctuation, the parallel derivative being taken 
along the total magnetic field comprising the background plus its perpendicular perturbation. 
Implicit in Hall-MHD is the assumption of cold ions, or at least r = Tj/Tg ^ 1. On the 
contrary, the gyrofiuid model is valid for arbitrary value of Tj/Tg. An other simplification 
can be reached in the limit k±pi ^ 1, i.e. ne = — 0/pf, in which case Eqs.(l)-(3) become 
similar to the reduced electron-MHD (see Appendix and Ref.[8]). 

Assuming that all fields vary like e*(^*~^-^), the linearized equations (l)-(3) for the Fourier 
components are 

Lorie — k\\kj_ip = 0, (9) 
u^^J + k\\{(f) - plrie) = 0, (10) 

He = (ro(6) - iWpI (11) 

Hence, this system yields the dispersion relation : 

with To (6) = e~^Io{b). The Fade approximation for the operator Fo(&) gives 

u;' = kl[l + klipl + p^)l (13) 

which is a standard expression for the frequency of the Alfven wave in the "kinetic" regime[l|. 

Let us notice that the gyrofiuid derivation of the the Alfven wave frequency gives the 
same result as its kinetic counterpart, however the latter also provides the imaginary part, 
i.e. the coefficient associated with Landau damping. In the model proposed below for 
the energy cascade produced by Alfvenic turbulence, any form of dissipation, including this 
collisionless one, is neglected at all except the smallest scales where it is supposed to balance 
the energy injected at large scales. For Alfvenic perturbation, relations between the density, 
the potential and the fiux function are provided by the equations (9)-(ll), the magnitude 
of the perpendicular magnetic field fiuctuation is given by B± = k^ip, the magnitude of 
the perpendicular component of the electric field is given by Ej_ = k±(j) and the parallel 
component by = p1k\\nf,. 

5 



III. ANISOTROPIC ALFVENIC TURBULENCE 



For Alfven waves, we observe from the energy integral Eq.(5), that the total energy density 
per wave-number in the perpendicular plane, i.e. the one dimensional energy spectrum 
denoted by E^^^, is 

Ek, = k^ij', (14) 

which coincides with the magnetic energy spectrum E^^ = 5Bj_/kj_. 

Locality of the non-linear interactions and constancy of the energy flux are assumed, 
hence the energy cascade rate e is 

6 = (15) 

tnl 

On the other hand, the non-linear time scale is given by 

rNL = ^ (16) 

with the generalized stream-function $ = (p — p^ng. Therefore, the non-linear time scale can 
also be written as a function of the magnetic perturbation, i.e. 

T^L = [Psklil + -^y^'k^^P]-\ (17) 
i — i 

with k±ij = {Ek^k±Y^^. From 1^, the energy spectrum is then obtained : 

T 

^ 

The Fade approximation for the latter expression gives 



E,,oce'/'pf/'kl'/\l + -^r/\ (18) 

i — i n 



E.,oc6^/^A://^[l + A:i(p^ + pD]-V3. (19) 

When Ti <^ Te, then Ekj_ oc e'^^^k_^_^^^{l + fc^p^)"^/^ showing a breakpoint for k±ps ~ 1. 
Hence, this is an Hall-MHD result. On the contrary, when <C Tj, then E^^ oc e'^^'^k]_^^^{l + 
k']_pf)^^^^ with breakpoint given by k±pi ~ 1. In general, the breakpoint occurs at k±p ~ 1 
with p = (p^ + PiY^"^- This breakpoint separates the MHD range with Ek^ oc e^/^A;^^/^ and 
the dispersive range with Ek^ oc e'^^^p^'^/^k_J^^ . 

Moreover, we have the following scaling relations: B± oc e^^^kj_^^^[l + k']_{p^ + pf)]'^^^^, 
Tie oceV3fci/'[l + fci(p2 + p2)]-2/3, OC e^/^ fc"'/' [1 + (p^ + p^)] "2/3 (1 + A^ip^) . Notice that 
He above refers to density fluctuations that arise from the kinetic Alfven wave compression. 
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The one dimensional energy spectrum for the perpendicular electric field is therefore 
El/k^ oc e2/3A;-^/^(l + fcip2)-4/3Q ^ klpff. For ~ T^, the latter relation is a good fit to 
some solar wind measurements |8|, [iSj. It corresponds to a power index —5/3 in the MHD 
regime when k±pi ^ 1 and —1/3 in the dispersive range when kj^pi ^ 1. When Tj ^ Te, 
i?_L oc e-'^/'^/c_|_^''^[l + k^plY^^ , which has only one break-point at kj^pi ~ 1. Notice however 
that in a plasma with Te ^ Ti then, Ei_ oc e-'^/^A;]^^^^[l + k\pf\'^'^/^{l + k\pj)^ which has two 
break-points at k^ps ~ 1 and at k^pi ~ 1. 

A theory for strong anisotropic Alfven wave turbulence was developed by Goldreich- 
Sidrarjl^. According to this theory, the anisotropy of the turbulence is fixed by the condi- 



tion, 



T, 



NLi 



(20) 



which is a balance between the linear and the non-linear dynamical time scales. This leads 
to the following scale dependent anisotropy relation : 



(21) 



This is also equivalent to the ordering relation 



5B^ ~ k\\/k^ < 1, 



(22) 



with 5-B_l(/c_l) given above. We end this section by emphasizing that the scaling relations 



for the energy spectrum and anisotropy in the dispersive scales 
of EMHD turbulence 16 



15| are similar to the ones 



IV. PARALLEL ELECTRIC FIELD SPECTRUM 



Alfvenic turbulence involves electric field fluctuations E which possess a component par- 
allel to the magnetic field B. We now discuss the spectral structure of the parallel electric 
field produced by Alfvenic turbulence. The magnitude of the parallel electric field is small 
compared to the perpendicular one, hence hardly measurable and not expected to change 
the total electric field spectrum. Nevertheless, the parallel electric field can be important 
for particle acceleration and cross-field transport induced by Alfvenic turbulence, or in fact, 
by any type of electromagnetic turbulence. The physical reason is that the parallel electric 
field can efficiently accelerate particles along the magnetic field lines which are bent in the 
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direction perpendicular to Bq due to the perturbation B^. Hence the correlation between 
acceleration and cross-field transport, and the pivotal role played by E\\. The results of the 
previous sections provide the scaling of the parallel electric field fiuctuation : 

E||«eV3p^fc||fcf [l + A;i(p^ + pD]-^/^ (23) 

together with Eg. (1211) for the relation k\\{k±). It is therefore easy to verify that the magnitude 
of ^11, function of /cy, scales like 

E\\{k\\) (X plkl (24) 

when <C e^/'^p"^/^, and like 

when k\\ ^ e^/'^p"^/^. Moreover, the magnitude of reaches its maximum. 

Eh - e^/^ /26) 

at the breakpoint k\\ ~ e^^^p^'^^^, which corresponds to to k±p ~ 1, with p = (p^ + pfY^"^-, 
where the two branches merge smoothly. It is clear that the condition kj^p ^ 1 does not 
have do be realized for E\\ to be of significant value because the latter reaches its maximum 
precisely at the boundary between the MHD and the "kinetic" scales. 



V. DISCUSSION AND CONCLUSION 



Based on a simple truncation of the gyrofiuid equations, we constructed a model for the 
energy cascade produced by Alfvenic turbulence which smoothly connect the large magne- 
tohydrodynamics (MHD) scales and the small "kinetic" scales. A similar model, with an 
emphasize on the effect of linear coUisionless dissipation was proposed in Ref.[20]. Scaling 
relations are obtained for the electromagnetic fluctuations as a function of k\_ and k\\ and 
particular attention is paid to the spectral structure of the parallel electric fleld produced 
by Alfvenic turbulence. The reason is that wave-particle interactions through the parallel 
electric force can produce anomalous transport in both velocity and coordinate space. A 
flrst principle understanding of such turbulent process is provided by quasi-linear theory 
applied to the gyro-averaged drift-kinetic equations. Within the framework of quasilinear 
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theory, the resulting diffusion equation for, say the electron distribution function / reads 
(see Ref.l) 

|-i:V|(;^)=f^(.-^„,,)4V/l. (27) 



with the operator V defined as 



V.„ + ^V.. (28) 



■'ce 



All quantities here are dimensional, f y is the component of the electron velocity parallel to 
the magnetic field and x is the coordinate transverse to Bq. The same type of equation can 
be written for the ion distribution function. A major ingredient entering the quasilinear 
diffusion coefficients, is the spectrum of the parallel electric field, which has been studied 
above for strong and anisotropic Alfvenic turbulence. The parallel electric field produced 
by the turbulence results in Landau damping of the turbulent field. The effect of any 
scale dependent damping process, can be accounted for on the energy cascade by writing a 
transport equation through wavenumber space as follow (see Ref. [20]): 

^ + ^'^l = -2^'^'- (29) 

at akx. tnl 

with 7(/c_L,fc||), representing the damping rate. Recently, the effect of Landau damping on 



the energy spectrum was studied this way, based on linear gyrokinetic calculations [20|] of the 
damping rate. However, the modification of the distribution function due to wave-particle 
interactions was not addressed. A possible approach is the standard quasi-linear theory. Let 
us notice that, in order for the fluid treatment to be valid, the number of accelerated particles 
must remain small compared to the number of core particles. Modeling of such a turbulent 
acceleration mechanism through the parallel electric field spectrum of Alfvenic fluctuations, 
its impact on the particles distribution function and back-reaction on the turbulent energy 
cascade, is the subject of ongoing work. 



VI. APPENDIX 



In this appendix we provide few additional details concerning the gyrofluid model. We 
also discuss some of its limiting cases (RMHD, ERMHD) and the corresponding model 
energy spectra. The normalized equations (l)-(2) follow from the electron continuity and 
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the parallel momentum equations, which can be written as 

dt- + ^[0, -] = d^v,, + V,,] (30) 

no -Do ^0 -Do 

5*^. + ^[^.A^] = -cd.^+^^d,^ + ^[^, A,] (31) 
Bo e no cBq uq 

with Vez = (c/47rno) V5_^z- This system is closed by the gyrokinetic Poisson equation 

ne = (ro-l)^0, (32) 

J- i 

Adopting the following "MHD" normalization, 

^ 7 ^ X , t X (j)C Az UeOOciTA. , ^ 

(t,x,0,A2,ne) = ( — ,-,- ), (33) 

ta L LvaBo LBo no 

with Ta = L/va, then Eqs.f l30|) - fl32l) become Eqs.(l)-(3). Hence, there are two non- 
dimensional parameters in this gyrofiuid model : say, ps and pi. When k±pi <C 1, the 
Poisson equation gives n^ = V^0, and therefore Eqs.(l)-(2) are equivalent to the following 
I0W-/3 Hall-MHD model : 

9fVl0+[0,Vl0]-V||J = O, (34) 
9tV^ + V||(p^Vi0-0) = 0. (35) 
In the MHD approximation, E\\ = and therefore the standard RMHD reads : 

atVi0 + [0,vi0]-V||j = o, (36) 

dtij - V||0 = 0, (37) 

Notice that the RMHD system is free of any non-dimensional parameter under the present 
normalization. The RMHD conserves the energy : 

E = J d^[{V±^y + {V±(j))% (38) 

Linearizing the RMHD system provides the frequency of the Alfven wave u = iA^y- More- 
over, since the relation = ±ip holds for Alfvenic fluctuations, it follows from the energy 
integral that the one-dimensional energy spectrum is 

Ek^ = k^iP\ (39) 
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which coincides with the magnetic energy spectrum, i.e. E^^ — SB^/k± since SB± — k±ip. 
Assuming locahty of the non-hnear interactions and constancy of the energy flux, the energy 
cascade rate e is 

6 = (40) 
tnl 

On the other hand, the non-hnear time scale is given by 

Hence, this time scale can also be written in term of the magnitude of the magnetic fluctu- 
ation, 

with k±tlj — {Ekj_k±y^'^. The expression for the energy spectrum in the MHD regime is then 
obtained : 

E,, oc e'/'kl'/'. (43) 

An other limit consists in taking k±pi ^ 1, in which case the Poisson equation (3) gives 
= —pfue, and hence Eqs.(l)-(2) become 

dtcl> + p'iViiJ = Q, (44) 

dt^ - (t^)V||(/> = 0. (45) 

An extension of this Electron-RMHD model, valid for a wider range of values of v^Jv\, was 
derived and studied previously in Ref.[8]. This ERMHD conserves the energy : 

/rfx[(V^^)^ + (^)S. (46) 

'Pi 

Linearizing the ERMHD system provides the frequency of the dispersive kinetic Alfven 



wave u) — ±k\\k±^T/l + rpi. Moreover, since the relation y^r/l -|- T{(f)/pi) — ±k±ip holds 
for kinetic Alfven waves, it is easily verified that the one-dimensional energy spectrum is 
again the magnetic energy spectrum Ekj_ = kjjp'^. The non-linear decorrelation time scale 
is given by 



1 + T 1 

(see Ref.[8]). This is indeed the k±pi » 1 limit of the relation (16), since Ue — —(t>/ pI also 
in this limit. Therefore, the expression for the energy spectrum in the dispersive range is 



E,, oc ^''pf'\^-^f''kl"\ (48) 

T 
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